LltlP Et = 









JLAB-THY-ll-1380 
NT-LBNL-1 1-013 
NT@UW-11-14 
UCB-NPAT- 11-008 
UNH-11-4 



SU{2) Low-Energy Constants from Mixed- Act ion Lattice QCD 

S.R. Beane/'2 W. Detmold,^'^ P.M. Junnarkar,^ T.C. Luu,^ K. Orginos,^'^ 
A. Parreno,^ M.J. Savage/ A. Torok,® and A. Walker-Loud^ 

^ (NPLQCD Collaboration) 

^ ^Albert Einstein Zentrum fiir Fundamentale Physik, 

^ Institut fiir Theoretische Physik, Sidlerstrasse 5, CH-3012 Bern, Switzerland 

^ ^Department of Physics, University of New Hampshire, Durham, NH 03824-3568, USA 

^Department of Physics, College of William and Mary, Williamsburg, VA 23187-8795. 
"^Jefferson Laboratory, 12000 Jefferson Avenue, Newport News, VA 23606. 

^N Division, Lawrence Livermore National Laboratory, Livermore, CA 94551. 

+^ ^Departament d'Estructura i Constituents de la 

1— H Materia and Institut de Ciencies del Cosmos, 

JlLi Universitat de Barcelona, E-08028 Barcelona, Catalunya. 

^ '^Department of Physics, University of Washington, Seattle, WA 98195-1560. 

I— I ^Department of Physics, Indiana University, Bloomington, IN 47405. 

I ^Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720. 

^ Abstract 

00 An analysis of tlie pion mass and pion decay constant is performed using mixed-action Lattice 

QCD calculations with domain- wall valence quarks on ensembles of rooted, staggered nj = 2 + 1 

MILC configurations. Calculations were performed at two lattice spacings of 6 ~ 0.125 fm and 

b ~ 0.09 fm, at two strange quark masses, multiple light quark masses, and a number of lattice 

volumes. The ratios of light quark to strange quark masses are in the range 0.1 < mi/mg < 0.6, 

• • while pion masses are in the range 235 < m.,^ < 680 MeV. A two-flavor chiral perturbation theory 

^ - — 

analysis of the Lattice QCD calculations constrains the Gasser-Leutwyler coefficients I3 and Z4 to 

^ be I3 = 4.04(40) (55) and I4 = 4.30(51) (gg). All systematic effects in the calculations are explored, 

d including those from the finite lattice space-time volume, the finite lattice spacing, and the finite 

fifth dimension in the domain-wall quark action. A consistency is demonstrated between a chiral 

perturbation theory analysis at fixed lattice spacing combined with a leading order continuum 

extrapolation, and the mixed-action chiral perturbation theory analysis which explicitly includes 

the leading order discretization effects. Chiral corrections to the pion decay constant are found 

to give /tt// = 1.062(26) (|g) where / is the decay constant in the chiral limit. The most recent 

scale setting by the MILC Collaboration yields a postdiction of = 128.2(3.6)(g;o)(3;3) MeV at 

the physical pion mass. 
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I. INTRODUCTION 



The masses and decay constants of the pseudo-Goldstone bosons are hadronic observables 
that Lattice QCD can now calculate with percent-level accuracy in the absence of isospin 
breaking and electromagnetism. This is primarily due to the fact that the signal-to-noise 
ratio of the ground state contribution to pion correlation functions does not degrade exponen- 
tially with time. While Lattice QCD calculations are still being carried out at unphysically 
large quark masses, with relatively coarse lattice spacings, and in modest volumes, chiral 
perturbation theory (xPT) can be used to describe the dependence of the pseudo-Goldstone 
boson masses and decay constants on these variables. Such a description involves a set of 
low-energy constants (LECs), which can be determined from experimental measurements, 
or from the Lattice QCD calculations themselves. The LECs that are extracted from the 
pseudo-Goldstone boson observables also appear in other physical processes, and therefore 
accurate Lattice QCD calculations of pion and kaon correlation functions arc beginning to 
translate into predictive power for other -more complicated- observables involving pions 
and kaons. 

xPT, the low energy effective field theory (EFT) of QCD, provides a systematic descrip- 
tion of low energy processes involving the pseudo-Goldstone bosons [1] . The theory consists 
of an infinite series of operators (and their coefficients, the LECs) whose forms are con- 
strained by the global symmetries of QCD. The quantitative relevance of these operators is 
dictated by an expansion in terms of the pion momentum and light quark masses suppressed 
by the chiral symmetry breaking scale, A^. At leading order (LO) in the two-flavor (n/ = 2) 
chiral expansion, the two coefficients that appear arc determined by the pion mass, m^r, and 
the pion decay constant, /j^. At ncxt-to-lcading order (NLO), there are four new operators 
in the isospin limit whose coefficients are not constrained by global symmetries [2]; these 
LECs are the Gasser-Leutwyler coefficients. Two of these LECs, h and Z2, can be reliably 
determined from low energy tttt scattering [3] . The LEG I3 governs the size of the NLO con- 
tributions to m^, while controls the size of the NLO contributions to f^^. Lattice QCD, the 
numerical solution of QCD, provides a way to constrain these coefficients, including those 
that depend upon the light quark masses. Further, as Lattice QCD calculations can be 
performed to arbitrary precision with appropriate computational resources, they will likely 
provide more precise determinations of the LECs than can be extracted from experimental 
data. A number of lattice collaborations have recently determined and I4 using Uf = 2, 
n/ = 2 + 1 and n/ = 2-|-l-|-l calculations of and with a variety of lattice discretiza- 
tions [4-11]. These efforts have been compiled into a review article [12] which has performed 
averages of these various computational efforts. It should be noted that there is an increas- 
ing number of Lattice QCD calculations performed at or near the physical point [6, 13-16], 
and it will be exciting to have reliable predictions of hadronic observables that do not rely 
on xPT. 

In this work, we focus on the determination of I3 and I4 from the pion mass and the pion 
decay constant using a mixed-action calculation with domain-wall valence quarks on gauge- 
field configurations generated with rooted, staggered sea-quarks. This serves to strengthen 
the case that the systematic effects arising from the finite lattice-spacing, which are unique 
to a given lattice discretization, can be systematically eliminated to produce results that are 
independent of the fermion and gauge lattice actions. There are already preliminary results 
from mixed-action calculations which can be found in Ref. [17]. 

Section II describes the details of the Lattice QCD calculation. In Sec. Ill, details of the 
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systematic uncertainties are presented. Continuum and chiral extrapolations of the results 
of the Lattice QCD calculations are detailed in Sec. IV. Conclusions are presented in Sec. V. 

II. DETAILS OF THE LATTICE CALCULATION AND NUMERICAL DATA 

The present work is part of a program of mixed-action lattice QCD calculations performed by 
the NPLQCD collaboration [18-32]. The strategy, initiated by the LHP Collaboration [33- 
38], is to compute domain- wall fcrmion [39-43] propagators generated on the Uf = 2 + 1 
asqtad-improved [44, 45] rooted, staggered sea quark configurations generated by the MILC 
Collaboration [46, 47], (with hypercubic-smeared [48-51] gauge links to improve the chiral 
symmetry properties of the domain- wall propagators). The predominant reason for the 
success of this program is the good chiral symmetry properties of the domain-wall action, 
which significantly suppresses chiral symmetry breaking from the staggered sea fcrmions 
and discretization effects [52-54]. This particular mixed-action approach has been used to 
perform a detailed study of the meson and baryon spectrum [37] including a comparison 
with predictions from the large- A^c hmit of QCD and SU{3) chiral symmetry [55, 56]. The 
static and charmed baryon spectrum were respectively determined in Refs. [57, 58]; the first 
caknilation of the hyperon axial charges was performed in Rcf. [59]; the first calculation of the 
strong isospin breaking contribution to the neutron-proton mass difference was calculated 
in Ref. [21], and the hyperon electromagnetic form factors were explored in Ref. [60]. The 
majority of calculations using this mixed-action strategy have been performed at only one 
lattice spacing, the coarse lattice spacing of 6 ~ 0.125 fm; a notable exception was the 
calculation of Bk [61], which included the fine MILC ensembles with b ~ 0.09 fm. In 
Ref. [62] , very nice agreement was found between the prediction of the scalar oq correlation 
function from mixed-action xPT (MAxPT) and the Lattice QCD calculations of the same 
correlation function [63]. This was an important check of the unitarity violations that are 
inherent in mixed-action calculations. 

A. The Lattice QCD Parameters 

In our previous works [18-32], on the b ~ 0.125 fm ensembles, domain- wall valence prop- 
agators were calculated on half the time extent of the MILC lattices by using a Dirichlet 
boimdary condition (BC) in the time direction. With the relatively high statistics that have 
now been accumulated, systematic effects from the light states reffecting off the Dirichlet 
wall are observed and are found to contaminate the correlation functions in the region of 
interest (see Fig. 1). This "lattice chopping" strategy has been discarded, and the valence 
propagators are now calculated with anti-periodic temporal EC's imposed at the end of the 
full time-extent of each configuration. The exception is on the heaviest light quark mass 
point of the b ~ 0.125 fm ensemble. At this heavy pion mass, the correlation function falls 
sufficiently rapidly to not be significantly impacted in the region of interest by the choice of 
BC. Further, this ensemble contributes very little to our analysis in Sec. IV. 

The parameters used in the present set of Lattice QCD calculations are presented in 
Table I. On the b ~ 0.125 fm configurations, light quark propagators computed by LHPC 
with anti-periodic temporal BC's are used for the three lightest ensembles [38]. Strange 
quark propagators are computed from the same source points in order to "match" the light 
quark propagators. In addition, calculations on the b ~ 0.125 fm ensembles with a lighter 
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TABLE I: The parameters used in the Lattice QCD calculations. 



6 ~ 0.125 fm ensembles 



3 


Ullbi 




L 


T 




L5 


6mf™^ 


6mp^ 


L dwf 

Ulllg 


J,™ res 


src ^ ^^cjg 


6.76 


0.007 


0.050 


20 


64 


1.7 


16 


0.0081 


0.001581(14)" 


0.081 


000895f3) 


4 X 468 


6.76 


0.007 


0.050 


24 


64 


1.7 


16 


0.0081 


0.00164(3) 


0.081 


0.00091(2) 


8 X 1081 


6.76 


0.010 


0.050 


20 


64 


1.7 


16 


0.0138 


0.001566(11)" 


0.081 


0.000913(2) 


4 X 656 


6.76 


0.010 


0.050 


28 


64 


1.7 


16 


0.0138 


0.001566(11)" 


0.081 


0.000913(2) 


4 X 274 


6.79 


0.020 


0.050 


20 


64 


1.7 


16 


0.0313 


0.001227(11)" 


0.081 


0.000836(3) 


4 X 486 


6.81 


0.030 


0.050 


20 


32 


1.7 


16 


0.0478 


0.001013(6) 


0.081 


0.000862(7) 


24 X 564 












6« 


0.09 fm ensembles 








/? 


bmf^ 


bmf^ 


L 


T 


Ms 


L5 


brnf"^ 








Nsrc X Ncfg 


7.08 


0.0031 


0.031 


40 


96 


1.5 


40 


0.0038 


0.000156(3) 


0.0423 


0.000073(2) 


1 X 170 


7.08 


0.0031 


0.031 


40 


96 


1.5 


12 


0.0035 


0.000428(3) 


0.0423 


0.000233(2) 


1 X 422 


7.09 


0.0062 


0.031 


28 


96 


1.5 


12 


0.0080 


0.000375(4) 


0.0423 


0.000230(3) 


7 X 1001 


7.11 


0.0124 


0.031 


28 


96 


1.5 


12 


0.0164 


0.000290(3) 


0.0423 


0.000204(2) 


8 X 513 



Provided by LHPC [38]. 



than physical strange quark mass have been performed. Statistics on three h ~ 0.09 fm 
ensembles have been accumulated, with the lightest pion mass being ~ 235 MeV. Finally, 
approximately 6500 trajectories have been completed on an additional rooted staggered 
ensemble with the MILC parameters 

13 = 6.76, 6mP = 0.007, 6mf^ = 0.050, ^ = 24^ x 64 , (1) 

and measurements have been performed on them. 



B. Results of the Lattice QCD Calculations 

Correlation functions with the quantum numbers of the 7r+ were constructed from propa- 
gators generated from a smeared source with both smeared (SS) and point (SP) sinks. To 
determine the pion mass, the correlation functions were fit with a single cosh toward the 
center of the time-direction. 

C(sx){t) ~ A^sx) e-"^-^/' cosh(m,(t - T/2)) , (2) 

where X = S,P. Fits incorporating excited states over larger time ranges produced consis- 
tent results for both and A(^sx) ■ The pion decay constant was then determined from the 
extracted overlap factors, A(^sx), along with the input quark masses and computed values 
of the pion mass and residual mass [64] , using the relation 

In the limit L5 ^ cxo, the residual chiral symmetry breaking in the domain- wall action 
vanishes and mp® — > 0. In addition to these valence quantities, the mixed valence-sea pion 
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FIG. 1: EMP's of the pion correlation functions on the b 0.125 fm ensembles. For comparative 
purposes, the effective masses from the correlation functions with Dirichlet EC's in time are shown 
for the lightest ensembles (slightly offset for visibility). 



correlation functions have been calculated to extract the mixed-meson masses, as was done 
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FIG. 2: EMP's of the pion correlation functions calculated on the large volume b ~ 0.125 fm 
ensembles. For comparative purposes, the effective masses obtained in the smaller volumes are 
shown (slightly offset in time for visibility). 



TABLE II: The pion masses and decay constants from the Lattice QCD 
uncertainty is statistical and the second is systematic determined from the 



calculations. The first 
fit range. 





X T^al X L5 


brun 


bfn 






m007m050 


20^ X 64 X 16 


0.18159(42)(g) 


0.09293(45)(|i) 


0.2553(15) 


3.63 


mOlOmOSO 


20^ X 64 X 16 


0.22298(26)(|^) 


0.09597(27)(^?) 


0.2842(15) 


4.46 


m020m050 


20^ X 64 X 16 


0.31091(27)(f[]) 


0.10204(26)(i?) 


0.3516(09) 


6.22 
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0.37469(22)(20) 
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40^ X 96 X 40 
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0.1344(14) 


4.13 


mOOSlmOSl 
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4.06 
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0.14530(15)(ii) 


0.06539(14)(i^) 


0.1632(10) 


4.07 


m0124m031 


28^ X 96 X 12 


0.20043(17)(}^) 


0.07032(19)(1[]) 


0.2153(03) 


5.61 



in Ref. [65]. 

The results of the Lattice QCD calculations are given in Table II. The quoted fitting sys- 
tematic uncertainties are determined by varying the fit range. Effective mass plots (EMP's) 
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FIG. 3: EMP's of the pion correlation functions on the b ~ 0.09 fm ensembles. 



for the full-volume correlation functions are generated with a cosh-style effective mass; 



m'f = - cosh" 



C{t + r) + C{t-r) 
"'^ T [ 2C{t) 

while the others were generated with a log-style effective mass; 



(4) 



<=-ln^^j. (5) 

In Figs. 1-3 the EMP's of the SS correlation functions and the extracted pion masses are 
presented using r = 3. 

In Fig. 1, the effective masses from calculations with anti-periodic BC's imposed on the 
valence quarks, as well as those from the Dirichlet temporal BC's, are shown. Correlation 
functions from propagators generated with a Dirichlet EC (located at t = 22 and i = — 10 in 
the figures) show a significantly different behavior from those generated with anti-periodic 
BC's. It is for this reason that we have abandoned the Dirichlet BC in the generation 
of valence quarks. However, it is only the lightest ensemble on which the extracted pion 
mass determined with the Dirichlet BC is statistically discrepant from that generated with 
anti-periodic BC's. 

Interestingly, the correlation functions generated with anti-periodic BC's are not free of 
their own systematics. The EMPs exhibit an oscillation with a period of approximately 
1 fm, which is not simply explained by either the staggered taste-pion mass splittings or 
by the mixed-meson mass splittings. In the top panel of Fig. 2, the oscillations are more 
pronounced (with higher statistics). Comparing the EMP's from the b ~ 0.09 fm and 
b ~ 0.125 fm ensembles, the oscillations are seen to become more pronounced for lighter 
quark masses. As the statistics are increased, the amplitude of the oscillation becomes more 
significant and increasing L5 does not appear to ameliorate these effects. At this point, it is 
not clear if the oscillations are an artifact of this particular mixed-action, or originate from 
the domain-wall valence propagators. Similar oscillations are observed for calculations with 
domain- wall valence propagators computed on dynamical domain- wall ensembles, as shown 
in Fig. 11 of Ref. [10] and Fig. 2 of Ref. [66]. For the present work, the masses and decay 
constants are determined with fits that encompass at least one full period of oscillation, with 
the fitting systematic estabhshed through variations of the fitting ranges. 



C. Scale Setting 

To extrapolate the calculated pion masses and decay constants and make predictions at the 
physical pion mass, the scale must be determined. The MILC collaboration has performed 
extensive scale setting analyses on their ensembles, and it is used to convert the calculated 
pion masses and decay constants into ri units (extrapolated to the physical values of the 
light quark masses).^ In Table III these values are listed for the ensembles used in this 
work [47]. The MILC Collaboration has determined ri = 0.318(7) fm using the bb meson 
spectrum and ri = 0.312(2) (§) fm using to set the scale [47]. The value of 

ri = 0.312(2) (3) fm, (6) 



The distance ri is the Sommer scale [67] defined from the heavy-quark potential at the separation, 
r?F(ri) = -1. 
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TABLE III: ri/h from MILC [47]. The values extrapolated to the physical light quark masses (right 
most column) are used to convert from lattice units to ri units. 



ensemble masses 




{bmi,bms,P) 


''i{bmf'\bmt\l3) 


m007m050 
m010m050 
m020m050 

m030m050 


6.76 
6.76 
6.79 
6.81 


2.635(3) 
2.618(3) 
2.644(3) 

2.650(4) 


2.739(3) 
2.739(3) 
2.821(3) 
2.877(4) 


m0031m031 
m0062m031 
m0124m031 


7.08 
7.09 
7.11 


3.695(4) 
3.699(3) 
3.712(4) 


3.755(4) 
3.789(3) 
3.858(4) 


TABLE IV: The pion masses (normalized to the 
units. The third uncertainty is the systematic from 


light-quark masses) and decay constants in r\ 
the conversion to ri units. 


ensemble masses 


V 


(rim^)2 




m007m050 
mOlOmOSO 
m020m050 
m030m050 


20^ X 64 X 16 
20^ X 64 X 16 
20^ X 64 X 16 
20^ X 32 X 16 


9.310(43)(i?)(ll) 
8.861(21)(i)(10) 
8.384(14)(i0)(10) 
8.275(10)(?9)(12) 


0.2526(12)(ii)(03) 
0.2634(08) (fl) (03) 
0.2877(07) (g^) (03) 
0.3092(04) (10) (05) 


m007m050 
m010m050 


24^ X 64 X 16 
28^ X 64 X 16 


9.318(23)(ii)(ll) 
8.846(16)(lt)(10) 


0.2550(08)(ig)(03) 
0.2640(ll)(ig)(03) 


m0031m031 
m0031m031 
m0062m031 
m0124m031 


40^ X 96 X 40 
40^ X 96 X 12 
28^ X 96 X 12 
28^ X 96 X 12 


I0.123(62)(70)(ll) 
9.942(57)(i4)(ll) 
9.551(20)(20)(08) 
9.285(16)(ii)(10) 


0.2331 (45) (i) (03) 
0.2318(34)(38)(03) 
0.2477(05)(}?)(02) 
0.2713(07) (?^) (03) 



is used in this work to convert to physical units. 
III. LATTICE SYSTEMATICS 

In order to make contact with experimental measurements, the lattice QCD results must be 
extrapolated to the continuum and to infinite volume, as well as to the physical values of the 
light quark masses. xPT is the natural tool to perform these extrapolations, a consequence 
of which is that the LEC's can be determined. 

A. Light Quark Mass and Volume Dependence 

Generally, the chiral expansion at NLO involves analytic terms, chiral logarithms and scale- 
dependent LEC's. However, the pcrturbativc expansion can be optimized by setting the 
renormalization scale to lattice-determined quantities which vary with the quark mass, lead- 
ing to modifications at next-to-next-to-leading order (NNLO). For instance, the SU (2) chiral 
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expansion of m^r and /tt can be expressed as [12, 18] 



A = / 1-ein 



^Phy 



^Phy 



-Ah 



where 



and 



k = log 



A? 



phy\2 ' 



(mr) 



(7) 
(8) 

(9) 



and Aj is an intrinsic scale that is not determined by chiral symmetry. Here and /tt 
denote lattice- measured quantities, / is the chiral-limit value of the pion decay constant, 
and B is proportional to the chiral condensate. The "phy" superscript indicates that the 
relevant quantity is evaluated with the physical values of the pion mass and decay constant, 
for which we use the central values 



ffy = 130.4 MeV 



and 



m 



= 139.6 MeV . 



(10) 



In addition to the light quark mass dependence, the finite volume corrections to the pion 
masses and decay constants can be simply determined in the p-regime, defined by rriT^L ^ 1. 
At NLO in the chiral expansion, the finite volume corrections are given by [68, 69] 

= 87r^AiX(^, m^L) (11) 



^{FV)_ 



2Bmq 

A(^^)^ = -167r2AzX(e, m,L) 



where 



k(n) 



rriTrL 



n=l 



(12) 



(13) 



and k{n) is the number of ways that the integer n can be formed as the sum of squares of 
three integers, n = Yl^=i "^1 with e Z. 

The light quark mass dependence of m^r and /^r are known at NNLO in two-fiavor 
XPT [70] . In the ^ expansion, in infinite volume, they are 



2Bm„ 



1 + -« W (A- 

2 [ \^p^y 



16 1 - 9- - 7 
y + ^/i2-^/3-/4-^ln(e^^) 



f ln(0-^4ee"^ + f^M (14) 



and 



k 
f 



5 



L — In 



^Phy 



^ + l/-2-5r4--ln(e'^^) 
12 6 ^ 2 ^ 



+ 2/4ee"^ + ffcF 



(15) 



where Z12 = 7Zi + 8/2- 
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TABLE V: Expansion parameters mi/ms, ^, ^uix, ^sea - ^, Csea - ^ and 



'^sea 


V 


mi/ms 




qMix 


Ssea s 


Ssea ~ s 


rUq 


m007m050 


20^ X 64 X 16 


0.14 


0.0491 


0.096 


0.114 


0.0032 


0.165 


mOlOmOSO 


20^ X 64 X 16 


0.20 


0.0681 


0.111 


0.108 


0.0010 


0.102 


m020m050 


20^ X 64 X 16 


0.40 


0.1177 


0.150 


0.093 


0.0001 


0.038 


m030m050 


20^ X 32 X 16 


0.60 


0.1540 


0.186 


0.084 


0.0026 


0.021 


m007m050 


24^ X 64 X 16 


0.14 


0.0489 


0.096 


0.114 


0.0032 


0.165 


m0i0m(J5(J 


2(8-'' X G i X IG 


0.20 


0.0G7i 


0.111 


0.108 


0.0010 


0.102 


m0031m031 


40^ X 96 X 40 


0.10 


0.0360 


0.058 


0.050 


0.0004 


0.039 


m0031m031 


40^ X 96 X 12 


0.10 


0.0365 


0.058 


0.050 


0.0004 


0.109 


m0062m031 


28''^ X 96 X 12 


0.20 


0.0629 


0.079 


0.045 


0.0019 


0.045 


m0124m031 


28^ X 96 X 12 


0.40 


0.1037 


0.119 


0.038 


0.0054 


0.017 



B. Mixed Action xPT 

The low-energy EFT for mixed-action Lattice QCD calculations is well understood [52- 
54, 62, 63, 65, 71-78]. At NLO in the MA expansion, including finite volume effects, the 
pion mass and decay constant are given by 

- \ (Isea - e) [1 + In (0] - ll"" (Csea ^ + (^) 

+ 87r^AiX(e, m,L) + - i)^diX{m,L) , (16) 



2Bmq 



7 = i-«M.i.>fe)+a 



2 



where 



- (CMix - i) In (e^'O - (^sea "0+^4 (^) 

- 167r^ AiXdwix, m^Mi.L) , (17) 

AdtI{mL) = y k{n) ( K^{^mL) + K^{^mL) - ^^^^"7^^ ) (18) 

For the present calculations, the extra expansion parameters of the theory are defined as 

?Mix — p ^sea — p ^^ea — p ^^^^ 

where ^Tr^^^g is the taste-5 staggered pion mass, 6^Ai is the mass splitting of the taste 
identity staggered pion and b'^A'^^^ is the mass splitting of the mixed valence-sea pion [73, 78] , 
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TABLE VI: Finite volume corrections to m-j^ and /^r at NLO in MAxPT, as given in Eq. (20) and 
Eq. (21). For a quantity Y in the table, SY[FV]/Y = {Y[FV] - Y)/Y. 









b ^ 0.125 fm ensemble 




Quantity 


m007m050 


mOlOmOSO 


m020m050 


m030m050 




L = 20 


L = 24 


L = 20 


L = 28 


L = 20 


L = 20 


MAxPT: 5m^[FV]/m^ 


1.6% 


0.6% 


0.6% 


0.1% 


0.1% 


0.0% 


XPT: 5m^[FV]/m^ 


0.2% 


0.1% 


0.1% 


0.0% 


0.0% 


0.0% 


MAxPT: dU[FV]/U 


-0.3% 


-0.1% 


-0.2% 


-0.0% 


-0.1% 


-0.0% 


xPT: 6U[FV]/U 


-1.4% 


-0.5% 


-0.6% 


-0.1% 


-0.1% 


-0.0% 








6 « 0.09 fm ensemble 






Quantity 


m0031m031 




m0062m031 


m0124m031 






L = 40 




L = 28 




L = 28 




MAxPT: dm^[FV]/m^ 


0.4% 




0.4% 




0.1% 




xPT: Sm^[FV]/m^ 


0.1% 




0.1% 




0.0% 




MAxPT: 5U[FV]/U 


-0.2% 




-0.6% 




-0.1% 




xPT: dU[FV]/U 


-0.6% 




-0.9% 




-0.2% 





determined in Refs. [62, 65] and this work. In Table V, the values of the parameters relevant 
for the calculations arc listed. 

In analogy with finite- volume xPT, the pion mass and pion decay constant in finite- 
volume MAxPT are related to their infinite volume values at NLO via the relations 



f 1 °° 



k{n) 



2^ 2 

n=l 



4e 



Ki{^m^L) 



+ (Csea - ( KoiVnm^L) + K2{y/nm^L) - 2 



Ki{y/nm^L) 



, (20) 



and 



1 - 4CMix Yl ^(^) 



n=l 



(21) 



In the case of /^r, the finite- volume effects in MAxPT are somewhat suppressed compared to 
those in xPT. This is because the contribution from the "average" valence-sea type virtual 
pion in a one-loop diagram is smaller than from a valence-valence pion due to its larger 
mass [65]. In contrast, the pion mass receives a one-loop contribution from a hairpin dia- 
gram [79], which has enhanced volume effects compared to a typical one-loop contribution. 
In Table VI, the FV contributions to and from Eq. (20) and Eq. (21) are presented. 
On the hghtest two coarse ensembles, the NLO volume contributions to from MAxPT 
are substantially larger than those from xPT. Further, due to the high precision of the 
Lattice QCD calculations, the finite-volume volume contributions are larger than the uncer- 
tainties on the m007m050 ensembles. This is in contrast to the results of the Lattice QCD 
calculations of m7r, which show little volume dependence. In Ref. [80], it was demonstrated 
that NNLO xPT could increase the finite- volume contributions by as much as ~ 50% of the 
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— NLOMAyPT : 

- - NLO xPT : 
• ' m010m050 

0.003 0.004 0.005 



Exp[-m^L]/(m^L)-''''2 Exp[-m^L]/(m^L)-''/2 



FIG. 4: NLO finite-volume contributions, and an estimate of their uncertainty, in and 
MAxPT compared with the results of the Lattice QCD calculations on the m007m050 and 
mOlOmOSO ensembles. The central values have been chosen to coincide for the larger volume en- 
sembles. 



NLO contribution. In the case of MA^PT, with hairpin diagrams having enhanced volume 
effects, the importance of the NNLO contributions are likely to be even greater than in xPT. 
As these NNLO effects have not yet been calculated, the MA^PT finite-volume contribu- 
tions are assigned a 30% systematic uncertainty when performing the analysis in Sec. IV. 
In Fig. 4, the NLO finite-volume contributions in xPT and in MA^PT for the m007m050 
and mOlOmOSO ensembles are compared with the results of the Lattice QCD calculations. 
The xPT band is given by the range Am^r = (1 + 0.5)AmJ^'^, while the MA^PT corrections 
are given by Am^ = (1 ± 0.3) Am^^-^^""", where the central values have been chosen to coin- 
cide for the larger volume ensembles. The MA^PT finite volume contributions appear not 
describe the observed volume dependence of m,r, indicating the likely importance of NNLO 
contributions. In the case of /tt, the volume contributions are in good agreement with the 
results of the Lattice QCD calculations. 



C. Strange Quark Mass Effects 

The strange quark masses that are used in the present calculations are not tuned to 
the physical value [81]; the physical staggered strange quark mass was determined to be 
= 0.0350(7) and bm^^^ = 0.0261(5) on the b ^ 0.125 fm and b ^ 0.09 fm ensembles 
respectively [47]. By matching SU{3) xPT onto the two flavor theory, the implicit strange 
quark mass dependence can be determined and absorbed into the NLO LECs [82]; 

h{m„ mf^) = r3(mf ^) + 6h{ms, mf^) , dUm,, mf ^) = In (^-^) 

— — — — \ I 771 \ 

Ums. mfy) = Umfy) + 6Um,, mf >') , 6Um,, mf ^) = - In ( j (22) 
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These lead to mild corrections to I3 and U on both the coarse and fine ensembles, 

SI (m mP^n = I ^ ^ ^"""^ = 

" ' ^ ^ -0.02, 6 0.09 fm, 6mf^ = 0.031 ; 

rj I phy^ / 0.09, h ^ 0.125 fm, hmf^ = 0.05 

D. Residual Chiral Symmetry Breaking Effects 

The domain-wall action has residual chiral symmetry breaking due to the finite extent of the 
fifth dimension, L5, resulting from the overlap of the chiral modes bound to opposite walls 
in the fifth-dimension. The quantity m^^^ is the leading manifestation of this residual chiral 
symmetry breaking, and the effective quark mass of the Lattice QCD calculation becomes 

ruq = mf"f + m]""^ , (24) 

capturing the dominant effects of the residual chiral symmetry breaking appearing at LO in 
the chiral Lagrangian. However, it is known that there are sub-leading effects. Defining the 
quark mass through Eq. (24) and taking the standard definition of rrf^^ as the ratio of two 
pion to vacuum matrix elements [64] 

^'°""«Ry' '''' 

where Jg^ and J5 are pseudoscalar densities made respectively from quarks in the middle and 

boundaries of the 5*^^ dimension, the quantity nf^^ = m^'^^{bmi,b) depends upon the input 
quark mass and the lattice spacing (sec Rcf. [10] for a discussion of these effects). Conse- 
quently, the chiral Lagrangian receives a simple modification at NLO [83-85] . Following the 
method of Ref . [86] , the modifications to the chiral Lagrangian at NLO are 

7 res I ;res 

SjCres = ^ ^ M r {2BmgJ: + 2Sm^E^) tr {2Bm'^'J: + 2Bm'^'J:^) 

;res 

+ ^tr (a^Sa'^S^) tr (2Bm''"'^S + 2Bm'^^'^St) . (26) 
8 

The corrections to 111^^ and /tt arising from these new terms are 

k — ir and ^ = e — r , (27) 



2Bmq 2 ruq f ruq 



with 



327r2 

;res ^ ^res ^ (28) 

li 

where 73 = —1/2 and 74 = 2 [2]. As with the coefficients l\, these l"^^^ coefficients are not 
universal and depend upon the choice of lattice action used. 
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TABLE VII: Parameters used to isolate nf^^ effects. The L5 = 16, 24 calculations were used to 
tune the quark mass for the L5 = 40 calculation in such a way that the sum h(mi + nij^^) was the 
same (within ~ 0.7%J for the L5 = 12 and 40 calculations. 



Ensemble 




hmi 






6m7r 


bfn 


4096f21b708m0031m031 


12 


0.0035 


0.000428(03) 


0.109(1) 


0.10160(22)(|^) 


0.0617(12)(iO) 




16 


0.0030 


0.000321(11) 


0.0987(3) 








24 


0.0030 


0.000229(12) 


0.071(4) 








40 


0.0038 


0.000156(03) 


0.039(1) 


0.10328(32)(36) 


0.0621(09)(i^) 



The new operators in Eq. (26) were found to give the dominant uncertainty in the pre- 
diction of the / = 2 vrvr scattering length at the physical pion mass [25] as the l^^^ were 
unknown. Therefore, for tttt scattering, and for other observables, it is important to deter- 
mine the 11^^, which can be done simply by performing calculations with different values of 
L5 on the same ensemble . The fine MILC ensembles, with b ~ 0.09 fm, at the lightest quark 
mass point were used to perform calculations with L5 = 12 and L5 = 40. The quark mass, 
defined by Eq. (24), was tuned to be the same for both Ls's, which was achieved to within 
0.7% accuracy (giving the same value of up to ~ 3%). The results of the calculations 
are presented in Table VII. The values of and ^4^^ that are determined by the Lattice 
QCD calculations are presented in Sec. IV. 

IV. CHIRAL, CONTINUUM AND VOLUME EXTRAPOLATIONS 

The chiral and continuum extrapolations of the results of the Lattice QCD calculations 
are performed in two different ways. The first method is to fit the LEC's of xPT to the 
b ~ 0.125 fm and b ~ 0.09 fm calculations independently. The extracted LECs are then 
extrapolated to the continuum limit, using the ansatz^ 

A(6) = Ao + A2 (^^^ . (29) 

This analysis is performed at both NLO and NNLO in the chiral expansion. The second 
method to perform the continuum and chiral extrapolations is to use MA^PT, which leads to 
determinations of the LEC's that are consistent with those obtained with the first method. 
This lends confidence that the discretization effects are small enough to be captured by the 
MA^FT formulation. 

Before proceeding, it should be noted that the light quark masses are given in lattice 
units and have not been converted to a continuum regularization scheme. As the product 
rrLgB is renormalization scheme and scale independent, the values of the LEC B, which we 
determine, have not been properly converted to a continuum regularization scheme. For this 
reason, we do not provide the results of this quantity. 



^ The leading discretization corrections in the current formulation of MA lattice QCD scale as OifP'). 
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TABLE VIII: Results of the fixed lattice spacing NLO analysis of m.,^. Max mi/ms denotes 
the maximum value of the ratio of light quark masses used to perform the analysis. 



TV It 

Max 




0.125 im 






mi/ms 


h 


^stat+sys 


dof 


Q 


0.4 


5.09(06)(52) 


18.1 


3 


0.00 


0.6 


4.60(03) (36) 


46.6 


4 


0.00 






b « 0.09 fm 






ml iris 


h 


^stat+sys 


dof 


Q 


0.4 


4.05(10)(40) 


3.31 


1 


0.07 



TABLE IX: Results of the fixed lattice spacing NLO X-PT analysis of f^^. Maxmi/mg denotes the 
maximum value of the ratio of light quark masses used to perform the analysis. 



Max 




b ~ 


0.125 fm 






mi/nis 


nf 


k 


^stat+sys 


dof 


Q 


0.4 


0.2166(10)(40) 


4.78(06) (20) 


2.35 


3 


0.50 


0.6 


0.2109(07)(13) 


5.28(03)(10) 


15.3 


4 


0.00 






b ~ 


0.09 fm 






mi/nis 


rif 


h 


Xstat+sys 


dof 


Q 


0.4 


0.1983(16)(34) 


5.48(13)(28) 


0.15 


1 


0.69 



A. Method 1: xPT and Continuum Extrapolation 
1. NLO SU{2) 

The pion masses and decays constants obtained in the Lattice QCD calculations on the 
b 0.125 fm and b ^ 0.09 fm ensembles are used to determine the LEC's at NLO in xPT by 
independently fitting to the expressions in Eq. (7) and Eq. (8), including the FV corrections 
in Eq. (11) and Eq. (12). Strange quark mass effects are included by using Eq. (22), but 
residual chiral symmetry breaking effects, such as those described by Eq. (27), are not. Both 
the mass and decay constant depend upon two LECs each, as seen from Eqs. (7) and (8) 
Including the larger volume calculations, the complete set of results presented in Table IV 
utilizes six data sets on the b ~ 0.125 fm ensembles and three on the b ~ 0.09 fm ensembles. 
For each of the NLO fixed lattice spacing fits that are presented in Table Vlll and Table IX, 
the maximum value of mi/ms used in the fit is listed. On the b ~ 0.125 fm ensembles, the 
ratio is in the range mi/mg — 0.14 — 0.6, while on the b ^ 0.09 fm ensembles the ratio is in 
the range mi/nis — 0.1 — 0.4^. 

From the quahty of fit given in Tables VIII and IX, it is clear that the NLO xPT formula 



In addition to giving the a-nd the number of degrees of freedom {dof) in the fit, the Q-value, or 
confidence of fit, is also provided, 

Q= r dx^r{x\d), (30) 
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TABLE X: Results of the continuum NNLO xPT analysis of uItt and /, 



Max 






b « 0.125 fm 










mi/nis 


rif 


^3 


U kM 


kp 


Xstat+sys 


dof 


Q 


0.4 


0.233(04) (08) 


7.95(35)(60) 


2.63(37)(67) 29(3)(4) 


21(6)(10) 


0.53 


4 


0.74 


0.6 


0.230(02) (03) 


5.83(14)(18) 


2.95(14)(24) 14(1)(1) 


16(2)(3) 


10.0 


6 


0.12 








b K 0.09 fm 












nf 




h kM 


kp 


Xstat+sys 


dof 


Q 


0.4 


0.203(11)(15) 


5.61(67)(73) 


4.1(1. 1)(1.6) 19(5)(5) 


2(17)(25) 











for m,r fails to describe the results of the Lattice QCD calculation at either lattice spacing, 
while the NLO xPT formula for f^^ describes the results on the lightest three h ^ 0.125 fm 
ensembles well and describes all the results on the h ^ 0.09 fm ensembles. Taking the 
results of the fits with mi/ms < 0.4, a continuum extrapolation of the extracted LECs using 
Eq. (29) gives 

[3 = 3.2(0.2)(1.2) and [4 = 6.3(0.3)(1.1) . (32) 

The NLO xPT determination of I3 must be taken with extreme caution (and essentially 
discarded) as the fit to is poor. This (relatively) large value of I4 extracted at NLO is 
consistent with the JLQCD NLO results using Uf — 2 overlap fermions [5] . 



2. NNLO SU{2) 

The pion mass and decay constant at NNLO in xPT, given in Eq. (14) and Eq. (15), depend 
upon two additional LECs, /cm and kp, in addition to the appearance of further NLO LEC's 
li2 — 111 + 8/2- Both Zi and I2 are reasonably well determined from tttt scattering [3], 

= -0.4(6) and I2 = 4.3(1) . (33) 

To perform the fits at NNLO, these values of h and I2 are used as input. Normal distributions 
of li and I2 are generated with means and variances given by Eq. (33), which are then 
used in the fitting process. This allows for a determination of the systematic uncertainty 
generated by their use as input parameters. In fitting to the results of the calculations on 
the b ~ 0.09 fm ensembles, there are six Lattice QCD results, and six fit parameters. The 
results of this analysis arc collected in Table X. The NNLO xPT is found to describe the 
results of the Lattice QCD calculations for both m^r and f„. Taking the b ~ 0.125 fm and 
b ~ 0.09 fm fit and using them to perform a continuum extrapolation, 

[3 = 3.3(1.4)(1.7) and [4 = 5.8(2.4)(3.5) (34) 



where 

^fe^*=5i7^(x=)'""e-«''^ (31) 

is the probabiUty distribution function for with d, degrees of freedom. (The Q-value represents the 
probabihty that if a random samphng of data were taken from the parent distribution, a larger would 
result.) 
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TABLE XI: Fit ranges used in the MAxPT analysis. For a given fit, A-E, the maximum value of 
mi/nis (sea-quark masses) is given. 



Fit Max mi/mg 

COARSE COARSE FINE 
L = 20 L = 24, 28 



A 


0.20 


0.20 


0.20 


B 


0.20 


0.20 


0.40 


C 


0.40 


0.20 


0.20 


D 


0.40 


0.20 


0.40 


E 


0.60 


0.20 


0.40 



TABLE XII: Results from NLO MAxPT fits to (rim^)^ / {nmg) . 









LECs 










Fit 


^3 


lb 


jres 




^stat+sys 


dof 


Q 


A 


4.27(23) (i) 


-1.23(21)(i) 


i4(6)a) 


-0.6(L6)(i|) 


L41 


2 


0.49 


B 


4.11(21)(i) 


-1.09(19)(i0) 


19(5)(i) 


-2.9(0.9)(2:0) 


2.33 


3 


0.51 


C 


4.10(19)(ii) 


-1.16(20)(iO) 


17(6) (i) 


-1.4(1.5)(f:f) 


1.78 


3 


0.62 


D 


4.10(19)(i) 


-1.09(19)(^9) 


19(5) (E) 


-2.8(0.8)(ii) 


2.33 


4 


0.67 


E 


4.10(19)(i) 


-1.13(18)(i«) 


18(5) (i) 


-2.7(0.7) (1;^) 


2.36 


5 


0.80 



are obtained, consistent with those from the NLO analysis. These results must also be 
treated with caution due to the small number of calculations performed on the b ~ 0.09 fm 
ensembles. 

B. Method 2: Mixed Action xPT 

As in the continuum case, the rriT^ and analyses with MAxPT are decoupled at NLO in 
the expansion, but the results of the Lattice QCD calculations at both lattice spacings can 
be fit simultaneously. This allows for several choices of fit ranges, which are denoted as A-E 
in Table XI. The maximum value of m^/m^ used in the fits from the b ~ 0.125 fm and b ~ 
0.09 fm ensembles are hsted in Table XI. As discussed in Sec. Ill B, the NLO MAxPT volume 
contributions are assigned a 30% uncertainty as an estimate of NNLO effects. This additional 
uncertainty is combined in quadrature with the other quoted systematic uncertainties. 

1. NLO Mixed Action xPT 

Fits are performed over the ranges listed in Table XI, the results of these analyses are 
collected in Table XII and Table XIII. There are a few observations to make. First, the 
NLO MAxPT formula is capable of describing the results of the Lattice QCD calculations of 
mjr, unlike the NLO xPT formula. Second, the MAxPT provides a slightly better description 
of the pion decay constant than of the pion mass. In both cases, the NLO formulae is capable 
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TABLE XIII: Results from NLO MAxPT fits to nf, 









LECs 












Fit 


rif 


h 








Astat+sys 


dof 


Q 


A 


0.1847(61)(i0) 


5.80(52)(^|) 


0.6(0.9)(1;?) 


-2(12)(li) 


-3.8(5.5)(?j) 


0.27 


2 


0.87 


B 


0.1860(20) (1?) 


5.73(42)(i) 


0.5(0.8)(0:i) 




-2.7(2.6)(it) 


0.28 


3 


0.96 


C 


0.1812(26)(||) 


6.03(40)(i) 


0.8(0.8)(?:8) 


-5(12)(}f) 


-6.1(4.4)(i:3) 


0.32 


3 


0.96 


D 


0.1841(17)(i) 


5.99(39)(1?) 


0.4(0.8)(0;i) 




-0.9(2.4)(|f) 


0.58 


4 


0.97 


E 


0.1797(12)(if) 


6.10(40)(i) 


0.9(0.8)(0;i) 


-5(ll)(}^) 


-2.9(2.4) (1:1) 


3.48 


5 


0.63 



of describing the results of the Lattice QCD calculations over the full range of quark masses. 

As the Q-value has a probabilistic interpretation, it is convenient to use it in forming 
weighted averages of the quantities that have been extracted with multiple fitting procedures 
and/or different numbers of degrees of freedom. For extractions of a parameter A from 
different procedures, each giving Aj with Qi, the weighted average 

A=^^, (35) 

can be formed. As each of the fits considered in this work, presented in Table XI, includes 
successively larger quark masses, this averaging will give more weight to the lighter quark 
mass values, where xPT is more reliable. Performing this Q- weighted averaging of the results 
from Tables XII and XIII gives 

UNLO] = 4.13(20)(^^) , k[NLO] = 6.09(40)(1^) , 

iriNLO] = 18(5)(^) , IriNLO] = -5(11)(}^) . (36) 

The value of ^3 is consistent with the average of all other Lattice QCD calculations [12]. 
However, the value of U is noticeably higher, but is consistent with that obtained with 
Nf = 2 overlap fermions and a NLO xPT analysis [5]. While the residual chiral symmetry 
breaking LECs are not well determined, they will help constrain the analysis of the / = 2 tttt 
scattering length [25] . 



2. NLO MAxPT + NNLO SU{2) xPT 

While the complete NNLO expressions for the pion mass and decay constant are not available 
in MAxPT, it is useful to consider the hybrid construction of NLO MAxPT plus NNLO 
xPT. As in the previous section, the NLO MAxPT volume contributions are assigned a 
30% uncertainty. Further, the infinite volume formulae for the NNLO contributions are 
used. While the fit values of the NNLO LECs will be polluted by discretization effects, the 
NLO Gasser-Leutwyler coefficients will be be free of these contaminations, and further, their 
extracted values should be stabihzed with the inclusion of these higher order contributions. 

The fit functions for and share two LECs. In principle, a correlated analysis 
should be performed, however, the correlations only exist at NNLO, and are expected to be 
insignificant. Results of uncorrelated fits are presented in Table XIV for the various data 
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TABLE XIV: Extracted values of the LEC's from NLO MAxPT plus NNLO xPT fitting of the 
Lattice QCD results. 









LECs 












Fit 


nf 


^3 


h 


ku 


kp 


Xstat+sys 


dof 


Q 


B 


0.186(9)(13) 


4.48(51)(f?) 


4.83(94)(|:|) 


13(5)(?) 


-8(17)(1I) 


2.22 


4 


0.69 


C 


0.188(7)(i?) 


4.12(30)(57) 


4.38(55)(i) 


8(2)(t) 


mill) 


2.17 


4 


0.70 


D 


0.193(5)(/o) 


4.00(28) (13^) 


4.10(44)(g) 


6(2)(i) 


5(6)GD 


2.99 


6 


0.81 


E 


0.194(3)( f) 


3.69(14)(i8) 


4.01(22)(36) 


3(1)(1) 


7(2)( I) 


3.63 


8 


0.89 




FIG. 5: I3 and I4 generated through a Monte-Carlo averaging of the fits in Table XIV. The 
histograms are generated with 10^ samplings. The vertical dashed lines represent the 16% and 84% 
quantiles. 



sets. Taking the Q-weighted average of these results gives 



/3[iViVL0]= 4.04(40)^5;, 
lT[NNL0] = n{b)O, 



h[NNLO]=A.'iO{^l){ll), 
ri'[NNLO] = 0(11) (12). 



(37) 



with /3[A^A^L0] and l^[NNLO] in good agreement with the averages given in Ref. [12]. At 
NNLO in the chiral expansion, corrections to the pion decay constant are found to be 



f 

j[NNLO] = 1.062(26)(^2 



40) 



(38) 



Setting the scale either by using r^^^ = 0.312(2) (I) fm from the MILC collaboration to 
determine f^^^, or by using the experimental value of f-,t+ to determine ri, gives 



/piiyfNNLO] = 128.2(3.6) (^:^)(y) MeV and rp[NNLO] = 0.306(9) (i^) fm. (39) 

where the last uncertainty in the postdicted value of /,r comes from MILC's determination 
of ri, Eq. (6). 

Figure 5 shows Monte-Carlo histograms of the extracted values of 1^ and U using the 
Q-weights to determine the ratio of samples to draw from each of fits A-E. The result of fit 
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FIG. 6: The result of NLO MAxPT plus NNLO xPT fit E described in the text, extrapolated to 
the infinite volume and continuum limits. The star denotes the experimentally determined value of 
f^+ (not used in the fitting). 



E for extrapolated to the infinite volume and continuum limits is displayed in Fig. 6. The 
inner (colored) band represents the 68% statistical confidence interval while the outer (gray) 
band results from the 68% statistical and systematic uncertainties combined in quadrature. 
The dashed vertical line is located at C,"^^^ determined from Eq. (10). 



C. Convergence of the SU{2) Chiral Expansion 

With the analyses performed in the previous section in hand, the convergence of the two- 
flavor chiral expansion can be explored. The resulting NLO and NNLO contributions to the 
quantities 

7-'' 

(both of which vanish in the chiral limit) are shown in Fig. 7. In both cases (the left and 
right panels of Fig. 7), it is the continuum limit and infinite- volume limit extrapolations 
that are displayed. In the case of m^, the NNLO contributions are negligible over most of 
the range of ^ used in our fits. Further, the total corrections to m^r are small, being less 
than ~ 15% over the full range of quark masses. In contrast, the corrections to /^r become 
substantial at the heavier pion masses, exceeding ~ 50% at the heaviest mass considered. 
Further, at the modest value of C,^ 0.08 the NNLO corrections become significant compared 
to the NLO corrections. 

In the left panel of Fig. 8, the determination of 1^ is shown. The results of the fixed lattice 
spacing xPT analysis from Sec. IV A 2 is displayed, as well as the continuum extrapolated 
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FIG. 7: The NLO and NNLO contributions to 



1 



J- (left panel) and — \ (right panel). 

Both of these quantities vanish in the chiral limit. The larger (red) dashed curves are the NLO 
contributions and the smaller (blue) dashed curves are the NNLO contributions. The solid (black) 
curve is the entire NLO + NNLO value. 
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FIG. 8: The present determination of I3 (left panel), and its comparison to the Lattice QCD 
average value [12] and phenomenological results (right panel). Some of the I3 results in the left 
panel have been given small offsets in for presentations reasons. 



value. Also shown are the values extracted from MAxPT at NLO, and from NLO MA^PT 
supplemented with continuum NNLO xPT, as discussed in Sec. IV Bl and Sec. IV B 2, 
respectively. The results of the MAxPT analyses are consistent with the continuum extrap- 
olated results, but with smaller uncertainties. This is not surprising as the mixed-action 
framework allows a simultaneous treatment of calculational results from multiple lattice 
spacings. This consistency lends confidence in the entire analysis. In the right panel of 
Fig. 8, the extraction is compared to the original estimates by Gasser and Leutwyler [2] as 
well as to the recent Lattice QCD average [12]. In Fig. 9, the analogous results for 1^ are 
displayed, although Ref. [12] does not provide an average value (citing insufficient reporting 
of the associated systematic uncertainties). 
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FIG. 9: The present determination of /4 (left panel), and its comparison with phenomenological 
results (right panel). (Ref. [12] does not currently provide a Lattice QCD average value for this 
quantity.) Some of the I4 results in the left panel have been given small offsets in (6/ri)^ for 
presentations reasons. CGL 2001 refers to Ref. [3]. 



V. RESULTS AND DISCUSSION 



We have performed precision calculations of the pion mass and the pion decay constant with 
mixed-action Lattice QCD. Calculations using domain-wall valence quarks and staggered 
sea-quarks were performed on a number of ensembles of MILC gauge-field configurations at 
different light-quark masses, different lattice-spacings, different volumes and different extents 
of the fifth dimension. The results of these calculations were extrapolated to the continuum, 
to infinite-volume and to the physical pion mass in two different ways, that were found to 
produce the same results within uncertainties. One method involved using two-flavor xPT 
to extract the LEC's, which implicitly include lattice-spacing artifacts. LEC's calculated 
at two different lattice-spacings were then extrapolated to the continuum. It is found that 
NLO xPT fails to describe the results of the Lattice calculations of m^, while NNLO xPT 
appears to be consistent with them. The second method was to use MA^PT where the 
lattice-spacing artifacts are explicit, and the extracted LEC's are those of the continuum, 
up to higher order contributions. A hybrid analysis was motivated to be sufficient, where 
the mixed-action NLO contributions were combined with continuum NNLO contributions 
to provide reliable extractions of the LEC's. These analyses have provided determinations 
of the Gasser-Leutwyler coefficients I3 and [4, 

[3 = 4.04(40) (1) and [4 = 4.30(51)0 (41) 

These values are consistent with the (lattice) averaged values reported in Ref. [12]. Our 
analysis also provides 

^ = 1.062(26)(1^), (42) 

which is to be compared to the lattice averaged value of /,r// = 1.073(15). Further, the 
extrapolated value of ri/^r and the experimentally measured value of /^+ provides a deter- 
mination of the physical scale ri, 

ri = 0.306(9) (J°) fm, (43) 



24 



which is to be compared with the MILC determination (on the same ensembles) of ri = 
0.312(2) (g) fm. It is interesting to note that, despite greatly enhanced statistics on the same 
ensembles of MILC gauge-field configurations, the uncertainty that we have obtained in the 
calculation of is somewhat larger than that obtained in Ref. [17]. 

The systematics in the calculations arising from the finite lattice volume and from resid- 
ual chiral symmetry breaking due to the finite fifth dimensional extent of the domain-wall 
action have been explored and quantified. Previously, residual chiral symmetry breaking 
contributions were identified to be the dominant source of uncertainty in Lattice QCD pre- 
dictions of the I — 2 nn scattering length [25] . While the present analysis has not been able 
to precisely determine these effects, the analysis resulted in constraints on the size of these 
contributions, 

ir^l7(5)(,l), Zr = 0(ll)(12), (44) 

which in turn can be used to reduce the uncertainties in the / = 2 tttt scattering length 
predictions. 

The predicted NLO mixed-action finite- volume contributions to the pion mass appear to 
be incompatible with the results of the Lattice QCD calculations, suggesting the importance 
of higher orders in the MA^PT expansion. A 30% systematic uncertainty is assigned to 
the NLO finite volume contributions to account for NNLO effects, leading to a consistent 
description of the results. 

In conclusion, we have found that a careful two-flavor low-energy effective field theory 
analysis of the Lattice QCD calculations of the pion mass and its decay constant can reliably 
determine the NLO Gasser-Leutwyler coefficients, [3 and J4, which are found to be in good 
agreement with the average of other determinations. In particular, mixed-action chiral 
perturbation theory which includes lattice-spacing artifacts explicitly, provides a reliable 
framework with which to perform chiral extrapolations of 111^^ and /^r to the physical light 
quark masses, and to determine /a and I4. 
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